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Abstract 

We study why it is quite so hard to make a superposition of superfiuid flows in a Bose-Einstein 
condensate. To do this we initially investigate the quantum states of N atoms trapped in a ID 
ring with a barrier at one position and a phase applied around it. We show how macroscopic 
superpositions can in principle be produced and investigate factors which affect the superposition. 
We then use the Bose-Hubbard model to study an array of Bose-Einstein condensates trapped in 
optical potentials and coupled to one another to form a ring. We derive analytic expressions for 
the quality of the superposition for this system, which agrees well with direct diagonalisation of the 
Hamiltonian for relatively small numbers of atoms. We show that for macroscopic superpositions 
to be realised there are essentially three straightforward requirements, other than an absence of 
decoherence, which become harder to achieve as the system size increases. Firstly, the energies 
of the two distinct superfiuid states must be sufficiently close. Secondly, coupling between the 
two states must be sufficiently strong, and thirdly, other states must be well separated from those 
participating in the superposition. 

PACS numbers: 03.75.Lm,03.75.Gg,03.75.Kk,03.75.Nt 
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I. INTRODUCTION 



One of the most puzzling questions in physics concerns why we do not see quantum su- 
perposition phenomena in the macro world [lj. In quantum mechanics superposition states 
are a crucial part of the description of the micro world, while in classical mechanics they 
are meaningless. This was highlighted by the thought experiment described by Schrodinger 
in 1935, where a cat is put into a superposition of alive and dead states [2j. Although 
Schrodinger's thought experiment is unachievable, there have been many more realistic pro- 



posals for creating systems that can be placec 
sufficiently decoupled from the environment 



in a macroscopic superposition if the system is 
3, 4]. In this paper we focus on what prevents 



the production of a cat state in superfluid flow, rather than on decoherence issues. Our 
discussion leads to some rather straightforward but nevertheless important conclusions. 
Superpositions have been experimentally studied in various forms [5J], including two slit 



electron interference and interferometry with atoms and molecules 

More recently larger systems have also been put into superposition states 



2] 



181 ] . Experimental signatures of large scale quantum superpositions have been observed in 



superconducting quantum interference devices (SQUIDS). The superposition was between 
states of different trapped flux or opposite currents flowing around the device. Truly macro- 
scopic systems can be defined as being visible with the naked eye or have a macroscopic 
measurable quantity associated with them. The currents measured in the SQUID consisted 
of approximately 10 9 Cooper pairs in a loop of diameter 140 /zm and produce a measurable 
magnetic flux, meaning; tunnelling between two macroscopically distinct states had been 



achieved 
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BECs hold promise for realising similar results as they contain 10 3 — 10 7 atoms with a high 
proportion in the same quantum state. They also have significant advantages over SQUIDs 
since they are highly controllable: the coupling between condensates and the strength of 
the interactions between atoms can be tuned over many orders of magnitude. They are also 
weakly interacting, which allows us to develop simp le models to investigate macroscopic 



quantum effects in BECs at the microscopic level [22j . To study macroscopic superpositions 
we need good control of the macroscopic parameters that influence the energy of the quan- 
tum states of the system. In the case of a SQUID the macroscopic variable is the external 
magnetic field, which can be applied accurately to the system over a large range of values. 
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This allows an easily controllable phase to be applied around the loop. By Larmor's theo- 
rem [23], the a»a,o gU e of app, y in g an external magn etic fie.d to charged pa.ie.es b a .oop 
is to move to the rest frame of the rotating potential that holds the neutral atoms of the 
BEC. The same effect could also be achieved by producing a flow of atoms round the loop 
by using Bragg scattering to imprint phases on the lattice sites 
effective magnetic field using two resonant laser beams 26, 2t| ] . 
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There have already been a number of theoretic al p roposals for creating cat states and 
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several using Bose- Einstein condensates (BECs) |28|, |29|, |30|, |3l|, |32j, |33j, [34]- In section [TT] 
we first look at a simple 1-D loop with iV particles on it and a barrier placed at one position. 
Using perturbation theory we show how a cat state of quasi-momentum (or flow) can be 
produced, and investigate the effects that inhibit its realisation. This analysis is developed 
in section III II using the Bose-Hubbard model (BHM) approach developed by Hallwood et 



al. 34J. It enables us to study quantum effects on an atomic level u sing the BHM 22| rather 
than the macroscopic approach used to describe the SQUID , l2lL |28| . The problem is 
tackled using the perturbation results derived in sect ion [IT] for an effective two state system. 
The original S x S matrix Hamiltonian, where S is the total number of states, is reduced to 
the desired 2x2 matrix leaving a single coupling term between the two distinct multiparticle 
states where the atoms share the same flow state. This is done by finding all the possible 
coupling paths between the two multiparticle states via intermediate states. From this model 
we are able to understand what factors inhibit the creation of a superposition as the system 
gets larger even if decoherence effects can be removed. In order to try and improve the 
"quality" of the superposition we consider differences between our model and that of the 
SQUID in section [TV] Our analysis supports that presented by Leggett in ref. [28( where he 
emphasises the role of the adiabatic limit in the operation of a SQUID. 



II. LOOP MODEL 



When discussing macroscopic systems we need to be careful with our terminology. In this 
paper we refer to a flow state as a state that describes the quasi-momentum or flow of each 
atom, *&kLk2,...,k N {rii ^2, • fiv), where and ki represent the position and quantised value of 
flow for atom i. The states obtained by diagonalising a Hamiltonian are called energy states 
with the ground state having the lowest energy. States with all atoms with the same flow 
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are called single flow states and the other states where not all atoms have the same flow are 
called multiple flow states. 
A cat state takes the form, 

*(r u r 2 ,...,r N ) = -^(* (n, r 2 , r N ) + e**i(n, r 2 , r N ))/y/2, (1) 
V2 

where \l/o( r i, r 2, ■■■, r N) and \l/ 1 (r 1 ,r 2 , ...,rjv) are both single flow states containing iV parti- 
cles. In these states all particles have the same flow, so only one flow is labelled. The most 
general description for a system of particles is the many-body Schrodinger equation, which 
has the form, 



1 N h 2 



2 ^ m. 



*fcifa...(r-i,r 2 ,...). (2) 



Here Viiti — rj) is the interatomic potential and V(ri) is the external potential. 

We cannot solve Eq. [2] directly, so we first simplify and then use perturbation theory 
to investigate some relevant features. Initially, we consider one particle confined in a one 
dimensional loop of radius R. There is only an angular dependence, 9, or if we wish to 
consider the position on the circumference, an x = R9 dependence. Eventually we will 
apply a phase around the loop, which can be considered as a rotation of the loop, then use 
perturbation theory to consider the effect of a barrier (see Fig. [T|). We then consider state 
with several particles. 

If all interactions and barriers are ignored the Schrodinger equation takes the simple free 
particle form: 

h 2 d 2 

ijj kn {x) = E kn ijj kn {x). (3) 



2m dx 2 

which has solutions ipk n {x) = ^e l2nknX ^ L , and energy E kn = j- (^) 2 k n 2 , where K = 

is a normalisation constant, L = 2-kR is the circumference of the loop and k n is an integer 

representing the different flow states. For example, when k n = there is zero flow, and 

when k n = 1 there is a 2tt phase around the loop indicating one quantum of clockwise flow. 

Throughout this paper we shall use the convention that positive phase variations correspond 

to clockwise flow. The quantisation of flow is due to the periodic boundary conditions of a 

loop. 

Applying a phase around the loop changes the momentum term in the Hamiltonian. This 
is similar to the vector potential term added to the Hamiltonian of an atom in an external 
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FIG. 1: Shows the system to be simulated. The particles are constrained to a loop of radius R and 
has a barrier centred at 4> = ir. 

I — I 

magnetic field or moving to the rest frame of a rotating object [23] . For instance if we rotated 
the loop and moved into its rest frame the momentum of an atom becomes p x — > p x — m\ x , 
where p x is the momentum of the atom in the lab frame and v x = Rd9/dt is the velocity 
that the loop is rotating in the lab frame. The additional velocity can be represented in 
terms of phase, 

v x = {h/m){<l>/L) (4) 
where <fi is the applied phase to the loop and L is the circumference of the loop. We can 



see where Eq. [4] comes from by considering the macroscopic wave function of a BEC 35] , 
ip(x) = ipoe l ®( x \ where is the phase of the condensate and ipo is the density, which we 
assume is constant for all x. When the phase is not constant throughout the condensate 
there is a velocity field associated with it, 

h f 
m J 

= ^ V$ (*)>- ( 5 ) 

If we move to the rest frame of a rotating loop then there will be an effective phase around 
the loop, 0, in addition to the phase due to the momentum of the atoms, By considering 
a linearly varying phase around the loop we get Eq. |H In the rotating frame the Schrodinger 
equation now takes the form: 

(Pa: - k(p/L) 2 

2^ ^ 0) = E kn i) kn 0) . (6) 

Even though the system is finite in size, an atom travelling around the loop will see a 
periodic potential, so the system is in effect a perfect infinite periodic potential. The solution 
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FIG. 2: Shows how the energy levels are affected by the phase around the loop. The dashed lines 
show a loop with no barriers, while the solid line shows a loop with a small barrier. C = S^- (^7) 2 
is a constant and a constant energy term has been subtracted from the solid lines so the graphs 
line up. 

of this equation is the Bloch function [36], ^k(x) = exp[iKx]uK(x), where K is the Bloch 
wave vector (or wave number in one-dimension). This produces the energy in the form, 




where in our model K = —<f)/L. The energy levels are shown in Fig. [2j The dashed lines 
show the calculated energy levels, where the zero flow state has zero energy at <fi = 0, the one 
flow state has zero energy at <fi = 2tt, etc... The solid lines show the energy levels when the 
ring has a barrier at some point. A constant energy term has been removed so the graphs 
line up. 

The solid lines were calculated using quasi-degenerate perturbation theory [3?! and show 
an anti-crossing in the energy levels at <p = it. In the system there are many flow states, 
but we will only consider the lowest two states that represent zero flow and one quantum 
of clockwise flow. We can ignore the other states if the coupling to them is small compared 
to their separation in energy. This is the case at <fi = it, where the ground state, to a good 
approximation, has the following form, 

^ = a il} (x) + ai^i(x) + corrections, (8) 

where ipo is an atom with zero flow, ip\ is an atom with one quantum of clockwise flow, and 
a and a\ are complex amplitudes that satisfy 0% + a\ ~ 1. In the absence of coupling the 
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FIG. 3: Shows a 2 level system where the energy levels have an anti-crossing at (f> = ir. E is the 
energy at the crossing, e is the energy of the unperturbed energy levels with zero at <j) = it, and 
Eq and El are the perturbed energy levels of the ground and excited states respectively. 

energy of the ipo(x) and ipi(x) states are given by Eq = E° + e(4>) and E® = E° — e{4>) where 
E° is the energy of both states at <fi = n and e describes the phase dependence of the two 
energy levels (see Fig. [3]). With coupling Voi the energies of the ground and first excited 
state become, 



The amplitudes ao and ai are obtained from diagonalising the Hamiltonian and a^ + a\ = 
1. The ratio of the amplitudes is [s7] , 



Sufficiently near = 7r we see e <C | Vqi | , so a /cii ~ 1 and we get a good superposition as 
long as Og + a\ ~ 1. e increases as we move away from (f) = n, so the zero flow or clockwise 
flow state's amplitude begins to dominate depending on whether the phase is increased or 
decreased. To make the superposition more stable we could increase the coupling between 
the single flow states, but this will increase the coupling to other states and dp + a? will no 
longer close to 1. The system can no longer be approximated to a 2- level system, and the 
simple form of the superposition will be lost. 

For two atoms, the full system, excluding spin, is described by the following Schrodinger 



E n = E°±(e(<f>) 2 + \V 01 \ 2 ) 1/2 , 



(9) 



Voi 



(10) 
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FIG. 4: Show how the energy levels are affected by the phase around the loop. C — h— ^7; 
constant. Some of the lines are now entirely degenerate with other lines. 



is a 



equation, 



Ek 1 ,k 2 ( ( / ) )' l f ) k 1 k 2 (Xl,X2) 



2m 2m 
+V I (x 1 -x 2 ) + V(x t ) + V{x 2 ) i'k 1 k 2 (xi,x 2 ), 



where Vj{xi — x 2 ) is the interatomic potential and V(xi) is due to a barrier. When Vi(x\ — 
x 2 ) = and V(xi) = 0, the wave function is separable and has the form, 



(12) 



with the associated energies, 



W« = ^( T ) 



(13) 



where N = 1/L. The energy spectrum is shown in Fig. HI For N atoms with the same flow, 
k, the energy is simply, 

This shows the gradient of the energy is proportional to the number of atoms in the system, 
N, and thus e in Eq. [TQl is also proportional to N. This means that as the number of 
particles increases, ao/ai will diverge from 1 more quickly as the phase is changed. For 
larger systems the range of values of <fi that still produce a cat like state will be narrow, 
which means it could become experimentally infeasible to produce a superposition. 
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For N > 1 there is no longer a simple crossing of two energy levels at = n. In the two 
atom case, the first degeneracy is between four flow states as shown in Fig. HI The solid 
lines are the single flow state i/)oo(xi,x 2 ) an d ^11 (^1 5 ^2), while the dashed line is both the 
multiple flow states (ipoi{ x i, ^2) ±^01(^2, x i))/y/2, which are degenerate for all <p. To create 
a superposition we need to increase the energy of the multiple flow states above the single 
flow states. This will reduce the amplitude of the multiple flow states in the ground state 
of the system. To do this we consider the interaction energy between the particles. 

To see how the desired separation might arise we assume the interaction energy is positive 
and it is the same between all particles. A system consisting of 2 identical bosons with 
different flow has wave function, 

i/tsih, k 2 ) = -= [V>(fci, k 2 ) + v>(fc 2 , ki)] , (15) 

where i[)(ki,k 2 ) — ^1(^1)^2(^2) and ij){k 2) k\) — '^(^OV'fciO 3 ^)- If the particles are in 
the same flow state then the wave function is simply ip(k, k). The expectation value for the 
interaction energy can be calculated using the interatomic potential V(|a;i— x 2 \). This gives, 

&(k,k)\V(\ Xl -x 2 \)\i;(k,k)) = V, (16) 

for all atoms in the same flow state. For atoms in different flow states we get, 

(Mh, k 2 )\V(\x 1 - x 2 \)\Mki, k 2 )) = k 2 )\V{\ Xl - x 2 \)\^(h, 

+&(k 2 , kx)\V(\xi - x 2 \)mk 2 , ki)) + k 2 )\V(\ Xl - x 2 \)\^(k 2 , h)) 

+ (^(k 2 , h)\V(\xi - x 2 \)mk u k 2 ))} = V + (y 12>21 + V 21 , 12 )/2, (17) 

where V = (if)(ki,kj)\V(\xi - x 2 \)\i>(k i ,k j )) and V ijM = (ip(ki,kj)\V(\xi - x 2 \)\tp(k k , ki)). 
The interaction energy has more terms when the atoms have different flow. If the interaction 
term is a short range 5-function, the energy shift is IV from the non-interacting case. We 
shall consider this case in more detail when we look at the Bose-Hubbard model. If we could 
use atoms with "medium" range interactions, e.g. dipolar ones 
additional flow state dependence. This is explored in section HVl 

As the number of atoms is increases more atoms can have different flow, which increases 
the energy shift from the non-interacting case. For three atoms the maximum energy shift 
is 6V from the non-interacting case, which happens when all atoms have different flow. If 
some of the atoms have the same flow the energy difference is smaller and the coupling to 
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the ground state is stronger. To produce a good cat state we need to separate the single flow 
states from the excited states. This is essentially what is realised in successful experiments 



to date. They achieve this in various ways. In the experiment described in Ref. it is 
the bound nature of the Cgo molecule that separates excited states from the two close lying 
states i.e. the molecule moving in two different directions. In the case of the SQUID the 
plasma frequency gives a gap to the excited states that have Cooper pairs with different 
flow [20] . The coupling to other flow states, which reduces the ability to make cat states, is 
similar to decoherence where the macroscopic state is coupled to other states via interactions 
with the environment. There is however a crucial difference, because the coupling is due to 
the system itself and can be reduced by making the barrier smaller, although this requires 
more accurate control of macroscopic parameters to make a cat state. 

From Eq. [10] we see we need |Voi| 3> e for ao/ai ~ 1 along with + a\ ~ 1. cIq + 
a\ decreases as the coupling to other states increases, because other flow states will have 
a significant amplitude. For N atoms with short range s-wave scattering the interaction 
potential has the form, 

N 

Vi = Vj2 5 (xi-Xj)- (18) 
The terms that affect the splitting are then given by, 



Voo-Vu = (Vo|VHVo) = 0, ( 19 ) 

L 

1 T O.TTX Ti -U -4-T »r 1 / T, C* / 

x 2 



v n r 

\V i\ = jjf-(N-l) J dx 1 ...dx N e 2 ^+- + *»» L 8(x 1 
N(N-l) 1 

= ^YT(2^' (20) 

where, tp = l/L N l 2 and fa = e^x/L^^/L ^J2-kx n /l j l n/2 _ Ag the n^ber of particles 
increases the (2tt) n - 1 term increases exponentially while the iV(iV — l)/2 term increases 
quadratically, so the overall coupling term, |Voi|, gets smaller. This makes it harder to get 
ao/ai ~ 1, therefore cat states become harder to make. 

The physical basis behind the effective coupling is the N(N — l)/2 ways of coupling 
from \N, 0} to |0, N). This coupling is, of course, effected directly by the energy defect 
to intermediate states. This makes the transition less likely for larger systems, because 
more energy is needed for all atoms to make the transition and gives the (2tt) n ~ 1 term, 
which decreases the energy gap at the anti-crossing. For example, to get from ^(0,0) to 
■0(1, 1) there is coupling through ip(0, 1), while the shortest path to couple from ip(0, 0, 0) 
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to 1, 1) is through ip(0, 0, 1) then ?/>(0, 1, 1). In the next section this will be investigated 
further with a more realistic system, the Bose-Hubbard model. 



III. COUPLED BECS IN A RING 



Hallwood et al. |34| presented a scheme for producing a multiparticle superposition of 
different superffuid flow states in a ring of coupled BECs using the BHM. In this scheme, 
it was possible to simulate the production of superpositions of all atoms stationary and all 
atoms rotating clockwise by applying a tc phase around a loop and changing the tunnelling 
strength between the sites. This was done using an adapted Bose-Hubbard model |35| that 



allows a change in applied phase between three sites that form a ring (see Fig. [5]). This gives 
the "twisted" Hamiltonian, 

H = - J[e^ /3 {a)b + tfc + Ja) + e"^ /3 (&t a + c% + o + c)] + U(a^a 2 + 6 1 " V + c^c 2 ), (21) 

where a, b and c are the annihilation operators of atoms on the three sites, U is the on-site 
interatomic interaction strength and J is the tunnelling strength between adjacent sites. The 



phase factors e ±4<?i//3 in the coupling terms are known as Peierls phase factors 35(. We note 
that does not have to obey the phase matching condition because it represents applied 
phase and not the phase of the condensate. 

Using this scheme we want to understand how the ground state is effected by experimental 
parameters and how cat states can be formed. The parameters we consider are the applied 
phase, 0, which needs to be ~ tt, the interaction strength, U, the coupling strength, J, and 
the number of atoms, N. The value of U / J determines the quantum phase of the system, 
giving a superffuid for U/J <C 1 and a insulator for U/J 1. We need to be in the near 
superfluid regime to achieve a good cat state, so atoms can tunnel freely, while there is still 
coupling between flow states. The state of the system can be found by direct diagonalisation, 
but there is a limit to how many atoms can be used of order 40 due to computational time. 
It is therefore necessary to find an approach that can look at higher numbers, which can be 
checked by comparison with exact calculations for lower numbers. 

It is convenient to consider a new complete orthogonal basis of operators that incorporates 
the phase matching condition, {a, (3, 7}. We will use the quasi-momentum basis, or flow 
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FIG. 5: The system consists of three sites, a, b and c, where atoms can tunnel between sites with 
tunnelling strengths J±, J2 and J3 and there is a clockwise phase difference of (ft/3 between each 
site. 

basis, because we will be producing superpositions of single flow states, 

a = (a + b + c)/v3, 

P = (a + be l2w/3 + ce U7r/3 )/V3, 

7 = (a + be~ l27T/3 + ce- l4lv/3 )/V3. (22) 

The new operator basis corresponds respectively to annihilation of an atom with zero flow, 
one quantum of clockwise flow and one quantum of anticlockwise flow. They obey the usual 
commutation relations and the quasi-momentum conservation rules for a periodic potential. 
Using these operators we can rewrite the "twisted" Hamiltonian in flow representation, 

H = -J{(2cv t a-/3 t /3-7 t 7 )cos(0/3) + \/3(^/3 - 7^) sin(0/3)} 
V + (3 ]2 (3 2 + 7 f V + 4(a t a/3 t /3 + 0+07+7 + 
+2(o 2 /3 t 7 t + /3 2 o t 7 t + 7 2 o+/3+ + h.c.}. (23) 

From Eq. fT23^) we see that the eigenstates of the system when £//J<Cl are just flow states. 
For U/J~ 1, there is coupling between the different flow states and when U/J^>1 the 
system is in the Mott regime where each site acquires the same number of atoms. 
We are interested in finding a superposition of the form, 

= a \N, 0,0) + a 1 |0,iV,0), (24) 
12 



where \N, 0, 0) is the state where all atoms are stationary and |0, N, 0) is the state where all 
atoms have one unit of clockwise flow. The terms in the ket represent the number of atoms 
in the a , j3 and 7 modes respectively. To study how good the produced cat state is, we 
calculate the ratio of the two amplitudes ao/ai as we did in section [III We consider systems 
where a 2 , + a 2 « 1, so a good cat state is formed when ao/ a i ~ 1- F° r a o + a i ~ 1 we require 
the system to be in the near superfluid regime, but we will see that this has a detrimental 
effect on the ratio ao/a\ « 1 when is not exactly tt. We require = tc for the energy of 
the flow states \N, 0, 0) and |0, N, 0) to be degenerate. It is, of course, important to see how 
small deviations from a 7r phase change will effect the quality of the superposition. This will 
enable us to assess the precision we shall need in an experiment. 

For large numbers of atoms the system states can be extremely complex. For now we 
will only consider atoms with all the same flow. In certain regimes this approximation is 
justified, because the atom-atom interactions raises the energy of multiple flow states faster 
than single flow states, so + a\ « 1. 

From Fig. [3] we see the flow state energy levels are given by E® = E° — e and E® = E° + e. 
Using Eq. EH we can reduce the problem to a 2 x 2 matrix Hamiltonian as in section [Til The 
energy and the amplitudes are described by Eq. [9] and Eq. [KB This shows the ratio of the 
amplitudes only depends on the coupling between the two flow states, Vqi, and the rate at 
which the energy levels change, e. For a stable superposition we would require Vqi ^> e, but 
for large V01 the perturbation must be big, so other states will begin to have a significant 
amplitude in the ground state. 

To use Eq. [10] we need to understand how e and <fi are related. The Hamiltonian for the 
case of no atom-atom interaction, i.e. U = 0, is diagonal in the flow basis, 

H = -J{(2a t a - ft (3 - 7 f 7 ) cos(0/3) - y/3(-^/3 + 7*7) sin(0/3)}. (25) 

The change in the Hamiltonian due to interactions is then given by, 

H' = -{a t2 a 2 + /3 t2 /? 2 + 7 t2 7 2 + 4(a + a/3 f /3 + ^7*7 + ^7) 
3 

+2(aW + /3 2 aV + 7 2 a t /5 t + h.c.}. (26) 

From this we can see that the unperturbed ground state energy has the form E\ = 
— 2 JN cos(0/3). If we subtract the energy of the two states at the crossing, = 7r/3, 
we get, 

e = E°(tt) - 2JiVcos(0/3) = JN - 2 JN cos(0/3). (27) 
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This shows that the gradient of the energy increases linearly with the number of atoms, and 
thus, the ratio ao/ai will diverge from 1 more quickly when moves away from ir as more 
particles are added. This will, of course, make it more and more difficult to produce a cat 
state. 

|0,0,4) 



1,0,3) |0,1,3) 




0,2,2) 



|0,3,1) 



3,0,0) |2,1,0) |1, 2, 0) |0, 3,0)|4,0,0) |3,1,0) |2,2,0) |1,3,0) |0, 4, 0) 



FIG. 6: Left: All the possible flow states for 3 atoms. The different colours show the coupling 
between different state. Right: All the possible flow states for 4 atoms. It is no longer possible to 
connect the pure flow states, so no superposition is possible between these states. 
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4,0,0) |3,1,0) |2,2,0) |1, 3,0) |0,4,0) 



FIG. 7: All the possible flow states for 4 atoms. The Hamiltonian that represents the coupling in 
this diagram has unequal tunnelling strengths, so coupling between the pure flow states is possible. 



From Eq. [23] we see there is no direct coupling from \N, 0, 0) to |0, N, 0), so we must 
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H 



(28) 



consider other flow states to calculate the total coupling through intermediate state. Note 
when U — there are no off diagonal terms and the coupling is not possible. In this case 
the eigenf unctions are just flow states. To understand what we mean by coupling path and 
intermediate states, lets consider a 4 state system where each state is represented by one of 
|1), |2), |3) and |4) with Hamiltonian, 

( ex V 12 Via V u \ 
V21 £2 V23 V24 
V31 V32 £3 V34 
\V 4 i V 42 V 43 £ 4 ) 

where e$ is the uncoupled energy of state i and Vjk is the coupling from state j to fc. Although 
there is direct coupling from |1) to |4) there are also indirect couplings through intermediate 
states 1 2) and |3). This is made more explicit in the four coupled equations that constitute 
the Schrodinger equation of the system. By eliminating the intermediate states we can 
simplify the problem to a two level system and a 2 x 2 matrix Hamiltonian. This gives a 
combined coupling from |1) to |4) of, 

T/ T/ V14V23V32 V12V24 V13V34 
^comb - ^14 ~ ~ ~ r + — r + 



(A-e 2 )(A-e 3 ) (A-e 2 ) (A-e 3 ) 



V12V23V34 V13V32V24 /^qN 

+ (\-e 2 )(X-s 3 ) + (A-e 2 )(A-e 3 )' 1 j 

where A is the ground state energy of the Hamiltonian that can be calculated using pertur- 
bation theory. If e 2 and £3 ^> e± and £4, and E\ ~ £4 then the ground state will have most of 
its amplitude in states |1) and |4), and a superposition between |1) and |4) is formed. What 
this shows us is the coupling is the addition of all the different paths from state |1) to |4) 
where the individual terms are given by the coupling to the intermediate states divided by 
(A-et). 

In the Bose-Hubbard model we can define a complete basis that describes all the different 
numbers of atoms with each flow where the total number of atoms is N. The coupling 
between these states is more clearly shown if we draw all the possible states in flow space 
and draw lines to show the possible coupling. This has been done for 3 and 4 atoms shown 
in Fig. [6] (a) and (b). In Fig. [6] (a) there is a coupling from |3, 0, 0) to |0, 3, 0) and |0, 0, 3) 
through |1, 1, 1), so under the correct conditions it will be possible to form a superposition 
between any two of the three single flow states. This happens when the phase is half 
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way between the phase of the two superposing states, so the energy of the two states is 
degenerate. In Fig. [6] (b) there is no path between the |4, 0,0), |0,4, 0) and |0, 0,4) states, 



so there is no coupling and, as has been shown [34j, no superposition is formed. This idea 
can be applied to higher numbers of atoms and it was found that coupling is only possible 
for commensurate numbers of atoms, i.e. the number of atoms divided by the number of 
sites must be an integer. 

This problem is resolved by adding different tunnelling strengths between the three sites, 
so the Hamiltonian now has terms coupling all adjacent flow states and Fig. [6] (b) becomes 
Fig. [7J This shows there are many different paths from one pure flow state to the other 
although the strength of the coupling depends on the differential tunnelling strength and 
atom-atom interaction strength. 

We were initially interested in seeing how this coupling changes as the number of atoms 
increases. Eq. [29] gave the form of the coupling from state |1), to |4) via states |2) and |3). 
Generalising to larger numbers of intermediate states we get, 

Vin = — . ,\ l ^ — . pn . + other possible paths. (30) 

(A - £j)(A - Sj)...(\ - Bp) 

All the different paths must be added together and each path has a combined coupling term, 
which is the product of the individual coupling terms along that path divided by (A — Sj) of 
each of the intermediate states, where Sj is the diagonal matrix element in the Hamiltonian 
of that intermediate state. As the number of atoms increases the length of the combined 
coupling increases and because is generally smaller than (A — Sj), so the overall coupling 
will decrease, even though the number of paths increases. From Eq. [10] we see that the ratio 
ao/ai will move away from 1 more rapidly for larger numbers of atoms as (f> deviates from 
7r. This will make it correspondingly harder to make a clean cat state. 

The combined coupling between the two single flow states has been calculated for iV = 3, 
6, 9 and 12 atoms and using Eq. [10] and [27] we have calculated a^/ai as a function of 
A(j) = <p — 7r as shown in Fig. [8] This can also be done for non commensurate numbers 
with unequal J, where there are many more coupling paths. Fig. [8] show the simulated and 
analytic results for the effect on a^/ai when the phase is not exactly ir around the loop. 
The analytic result for low numbers fit the simulated results very well. As the number of 
atoms increase the quality of the fit is reduced. This is due to the large number of states, 
so the total amplitude in the \N, 0, 0) and |0, N, 0) states is reduced. One can see in Fig. [8] 
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FIG. 8: Simulated and analytic results for the effect of changing the phase by A<p away from ir/3 
for (a) 3 atoms, (b) 6 atoms, (c) 9 atoms and (d) 12 atoms, where U/J = 0.1. 



the slopes drops more quickly as the number of atoms in the system is increased. This 
shows how experimentally difficult it will be to create a cat state, as the number of atoms 
is increased. 



IV. MODIFICATIONS TO THE MODEL 

The discussion in Sections [Til and II I II have demonstrated clearly that as the number of 
atoms is raised it becomes increasingly difficult to make a superposition. The analysis 
suggests that even for relatively moderate numbers it would be experimentally impossible 



to make superpositions. We know, however, that SQUID experiments 2CJ, |21[ have shown it 
is possible, so what feature is missing from our model? We believe that the crucial difference 
is our use of short-ranged interactions. In this section we shall, therefore consider the effect 
of long range interactions. We shall also consider briefly how having a continuous loop with 
a barrier, and not three sites, will improve our chances of making a cat state. 

We shall suppose that the long range interaction is due to either magnetic or electric 
dipole-dipole interactions between the atoms. An electric dipole-dipole interaction can be 
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induced by applying a strong d.c. electric field across the condensate [39(. It has also 
been shown to be possible to trap atoms with larger magnetic moments, so interactions 
beyond s-waved scattering must be considered 40]. There has in fact been a good deal of 



41 



42 



experimental and theoretical work on the dipole-dipole interaction 
this paper we shall use the formalism presented in the paper by Goral et al. 
magnetic dipole-dipole interaction potential is thus given by [411 ] . 

Ho /ii(r) ■ /i 2 (r) - 3(>i(r) • u)(/i 2 (r) ■ u) 
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45|. In 



dd 



„/ 1 3 



where the 



(31) 



An \r — r' : 

where u is the vector joining the interacting particles, fio is the magnetic permeability of 
the vacuum and fii are the magnetic moments of the atoms. This can be simplified by 
assuming the magnetic moments of each atom points in the same direction. Addin g th e 
s-wave scattering length potential, the total interaction potential then takes the form [38], 

3 cos 2 9 



Vi 



int 



Aixh 2 a,., A , 2 1 
-6{t - r') + d 2 - 



»/|3 



(32) 



m |r — r \ 

The first term is the usual short range interaction with an s-wave scattering length of a and 
atomic mass m. The second term is due to the dipole-dipole interaction of strength d. As we 
are only dealing with a 2-D plane we can set 9 = tt/2. The dipole-dipole interaction produces 
off site interactions in the BHM, so for the three site case there are nearest neighbour 
interactions only, and the Bose-Hubbard Hamiltonian becomes, 



H = -J[e^ /3 (a)b + tic + c f a) + e"^ /3 (tfa + c f 6 + a f c)] 

+U (a^ 2 a 2 + tf 2 b 2 + C t 2 c 2 ) + fA((at)V + {tffc 2 + (c*) V), 



(33) 



where a, b and c are the usual annihilation operators for the three sites, J is the tunnelling 
strength and Uq is the on-site interaction strength. The extra term U\ is the strength of the 
interaction between atoms on the adjacent sites. The interaction terms are given by, 



U„ 



Hr-r^lV^r-rOHr' 



iVrdV, 



(34) 



where |r,- — r~ 



47rcr/|k|, and i and j are the indices of the sites. w(r — r^) are Wannier 
functions used to describe the field operators in the Hamiltonian and are well localised on 
the lattice sites j^ . The two interaction terms are therefore given by, 

4Trh 2 a 



U 



rn 



\w(r - r )| 4 rf 3 r + d 2 
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|w(r-r )| 2 


w\ 


r'-r )| 2 




r — 


r' 


3 



d 6 rd i r 



3^' 



(35) 



and 

U l= <ff I^^-^IW-^IV ^. (36 ) 

The integrand in the second integral in Eq. [35] becomes infinite for r' = r, which is the 
case for on site interactions. This problem arises due to the pseudo-potential form of the 
dipole-dipole interaction. A better representation of the potential is required to solve this, 
which is investigated by Ronen et al. [46]. Even if the s-wave scattering is small and the 
dipole-dipole interaction is large we can see that U\ is always smaller than Uq. 

Again we are interested in the flow of the atoms, rather than position, so we transform 
to the flow basis defined by Eq. [22] Eq. [33] then takes the form: 

H = -J{(2a t a-/3 t /3- 7 t 7 )cos(0/3) + v^^-T^sin^/S)} 



(U + U x ) 



{a t2 a 2 + /3 t2 /3 2 + 7 t2 7 2 } + ^ t ^ {(ptafip + c*W7 + /Wt)} 



6 ' ' 111 6 

'^-{(aW + W + ^ + k.}. (37) 



6 

In order to get a good cat we need the extra long range interactions to increase the difference 
in energy between the single flow states and the multiple flow states. Examining the terms 
in Eq. [37] we can see the long range interactions reduce the coupling between the different 
states, but do not increase the energy difference between the pure and multiple flow states. 

Another difference between this model and that of the SQUID is the setup of the loop. 
Although our system forms a loop there are only three points at which the atoms can sit. 
In the SQUID there are sections that are superconducting. In these sections the long range 
interactions prevent atoms having different flow |47j, and make the cat state less phase 
sensitive. The gap to excited states in the charged superconductor is, of course, the plasma 
frequency. 

It is the junction that acts as a mechanism that creates the superposition of phases of 
the superfluid wavefunction. This is equivalent to the decaying atom in Schrodinger's cat 
thought experiment. Due to the small numbers of atoms in the junction it is relatively 
easy to create a superposition of states of uncertain phase, along with the demands of 
the uncertainty relation A<pAN m 1. The phase of the junction is coupled to the rest of 
the superfluid atoms, because of phase matching around the loop. If the loop contained a 
macroscopic number of atoms then this would in term form a macroscopic superposition. 
This is equivalent to the life of the cat being coupled to the atom. 
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In a SQUID, if the junction forms a superposition of phases then the rest of the loop will 
"follow" , as excitations of the super conductor are limited by the coulomb interaction, which 
produces a gap equal to the plasma frequency. Future work will examine this issue further 
to see how a SQUID like systems can be made using neutral atoms. Hopefully this will 
allow us to understand how to make stable cat states by showing the junctions only require 
numbers of atoms in the range we have investigated here and that they can be coupled to 
large numbers of superfluid atoms. Excitations are certainly a problem when the junction is 
coupled to the neutral superfluid regions. If the superposition is not produced adiabatically 
then phase and density fluctuations will be produced in the loop. Hopefully atom-atom 
interactions can be used to limit this effect in the same way that the coulomb force limits 
excitations in the SQUID 2^]. 



V. CONCLUSION 



We show three reasons, other than decoherence, why macroscopic quantum effects are 
so hard to see in Bose-Einstein condensates. We illustrate this by finding macroscopic 
superposition states of flow, firstly in a ID loop, then in a ring of coupled Bose-Einstein 
condensates. Firstly, we show that the energy of macroscopic states varies rapidly with the 
externally applied phase making it hard to keep the degeneracy needed for good coupling 
between them. Secondly, this coupling decreases as the size of the system increases, due to a 
larger number of intermediate states being involved. Finally, as more particles are introduced 
there are more possible flow states for the system to occupy. Even if the amplitudes for being 
in a multiple flow state are small it still reduces the total amplitude in the two macroscopic 
states. To avoid this we need a large energy gap between our macroscopic superposition 
states and the multiple flow states. Our results clearly show why it is hard to make systems 
with macroscopic quantum properties and will be useful in finding the best setup to produce 
a macroscopic superposition. These problem can be alleviated by coupling a microscopic 
mechanism that creates the superposition to a macroscopic system, as in Schrodinger's 
initial thought experiment where the microscopic mechanism is the decaying atom and the 
macroscopic object is the cat. In the context of the systems studied in this paper, an 
appropriate solution may be a microscopic Josephson junction coupled to a macroscopic 
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superfluid ring. 
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